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Abstract 

Two topics in heavy quark expansion are discussed. The heavy quark potential in 
perturbation theory is reviewed in connection to the problem of the heavy quark 
mass. The nontrivial reason behind the failure of the "potential subtracted" mass 
in higher orders is elucidated. The heavy quark sum rules are the second subject. 
The physics behind the new exact sum rules is described and a simple quantum 
mechanical derivation is given. The question of saturation of sum rules is discussed. 
A comment on the nonstandard possibility which would affect analysis of BR s i(5) 
vs. n c is made. 
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1 Introduction 



Heavy quark symmetry and the heavy quark expansion have played an important 
role in understanding weak decays of heavy flavors. Recent years witnessed signifi- 
cant success in quantifying strong nonperturbative dynamics in a number of prac- 
tically important problems via application of Wilson Operator Product Expansion 
(OPE). The advances were in theoretical understanding of strong dynamics of heavy 
quarks, as well as in practical applications to phenomenology of beauty decays. The 
discussion of both aspects can be found in the recent review M. It also included 
updates on such questions of general interest as the current theoretical uncertain- 
ties in extracting the KM parameters \V c b\ and \V u b\, our knowledge of the heavy 
quark masses and other basic parameters of the heavy quark theory, aspects of local 
quark-hadron duality relevant to the quest for \V c f,\ and \ V u b\. The interested reader 
can find there references to the original publications. 

In this contribution I address a few selected topics which have not been spelled 
in much detail in the literature. Even though some are of more theoretical nature, 

1 try to place emphasis on the general physical features understandable without 
presenting much of the technicalities. I discuss in some detail the notion of the 
heavy quark potential in QCD and its limitations, in particular as viewed through 
its relation to the heavy quark mass. The origin of the subtleties is explained which 
led to a failure in defining the low-scale running heavy quark mass in higher orders 
using the usual interquark potential. 

Another topic is the heavy quark sum rules, where the significant part is dedi- 
cated to new exact spin sum rules; applications to phenomenology are discussed. I 
conclude with the brief comment on charm counting in B decays. 

2 Heavy quark potential and the heavy quark mass 

Strictly speaking, the heavy quark potential appears in a somewhat different situa- 
tion compared to usual B decays, namely where heavy quark and antiquark (with 
small relative velocities) are present simultaneously. It attracted recently renewed 
attention in connection with the pair production of ti and bb near threshold, and 
through its connection with the problem of heavy quark mass. The review of appli- 
cations for the QQ system can be found, for example, in Ref. @. Here we give a 
more pedagogical discussion of the underlying problems. 

A closer look at the notion of the heavy quark potential in QCD reveals certain 
subtleties. It turns out that a literal analogue of potential interaction does not exist 
in QCD. 

The original notion of the potential refers to the interaction of infinitely heavy 
(static), or completely nonrelativistic heavy particles, which is instantaneous. The 
most familiar example is the electromagnetic interaction of heavy charges. The 
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Hamiltonian of such a potential system is given by 

i 1 

where are masses of particles and the potential V is a function of their instant 
coordinates. Taking the limit rrii — > oo (at fixed f{, which corresponds to semi- 
classically high excitation numbers of the quantum system in Eq. ([!])) eliminates 
quantum uncertainties in the coordinates and allows to measure the potential di- 
rectly as the position-dependent energy of the infinitely slowly moving collection of 
particles. This is well-known for QED where the potential of the charges q\ and q2 
(in units of the electron charge) is given by 

V^ D (R) = a em (0) ^ . (2) 

This expression is exact in the absence of light charged particles; known quantum 
corrections appear only if other matter fields are not much heavier than the scale 
l/R. 

The definition of the similar quantity in QCD turns out more tricky, since color 
sources are gauge-dependent. The color of the individual heavy quark remains fully 
quantum in nature and changes through interaction with gluons. In contrast to 
usual coordinates, the limit vtlq — > oo does not make color variables describing the 
state of the heavy quark semiclassical. 

To avoid this problem, the heavy quark potential is usually defined via the vac- 
uum expectation value of stretched Wilson loops, 

V(R) = - lim - In (Tr V exp [i I A^dxJ ) , (3) 

T^oo 1 \ Jc{R,T) ) 

where the rectangular contour C(R,T) spans distances R and T in the space and 
time directions, respectively. (It is usually assumed to be in Euclidean space as in 
Eq. (0)-) The gauge field is taken in the color representation of the heavy quark 
(fundamental representation for actual quarks). This definition is intuitively clear, 
since such Wilson lines describe propagation of infinitely heavy quarks. Moreover, 
such Wilson loops are readily computed in the U(l) gauge theory (free QED) and, 
of course, reproduce Eq. (|2]) (with qi = — ^2)- In QCD the perturbative expansion of 
V(R) has been computed through order Usually the potential in the momen- 

tum representation is considered, 

Vtf) -/** ViU) e * - ^g) 2 + ..] . (4, 

Taking the definition (^) for V(R), one faces a number of questions. What 
type of processes does V(R) incorporate and what are its properties, say, in the 
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perturbative expansion? Can it be used, say, in the Schrodinger equation similar to 
usual quantum mechanics or QED? It turns out that both questions are nontrivial 
and interrelated. 

There are important peculiarities in the thus defined static interaction in the QQ 
system; they were first analyzed by Appelquist et al. already in the late 1970's [|J. 
Due to the gluon self-interaction, the potential contains, in higher orders, the H- 
type diagrams shown in Figs. 1. Here the dashed line denotes the Coulomb quanta 
(they mediate instantaneous interaction in the physical Coulomb gauge), while the 
wavy lines are used for the transverse gluons. These diagrams propagate in time the 
transverse gluons. This means that at this level the problem ceases to be a two-body 
one, and includes more full-fledged dynamical degrees of freedom. 
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Figure 1: Examples of diagrams for heavy quark potential in QCD. Dashed lines are 
instantaneous Coulomb exchanges, transverse gluons propagating in time are shown 
by wavy lines, a) The convergent H-diagrams with the transverse gluon as a rung, 
b) Adding more Coulomb exchanges inside the "H" leads to infrared divergence in 
perturbation theory. 

Moreover, while the first H-diagram Fig. la appearing in order safely con- 
verges at the rung gluon momenta ~ 1/-R, including additional Coulomb exchanges, 
as in Fig. lb leads to an infrared divergence. With one exchange it is logarithmic; 
the formal degree of the infrared divergence increases with adding extra Coulomb 
quanta between the emission and absorption of the transverse gluon. Physics behind 
this infrared behavior was discussed in Ref . {4| : emission of the soft transverse gluon 
changes the overall color of the QQ pair and, therefore, modifies the interaction 
energy between them. The energy shift associated with the exchange of the trans- 
verse gluon depends nonanalytically on the energy denominator (in the language of 
noncovariant time-ordered perturbation theory), since the gluon can be arbitrarily 
soft. Formally expanding the exact result in a s includes expanding in this change 
of the Coulomb energy proportional to a s /R. Therefore, one obtains increasing in- 
frared singularities. These arguments suggested that resummation of the Coulomb 
exchanges in Fig. lb would render these diagrams finite. However, the effective in- 
frared cut-off is of the order of a s /R and, thus, the potential is not perturbatively 
infrared finite containing terms ~ 1/R ■ In a s starting at the order O [a s ) . 

This purely perturbative analysis shows that there is no direct analogue of the 
potential between heavy quarks in QCD. The QQ system incorporating all gluon 
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interactions is not a two-body system but includes actual propagation of gluons 
with energy small compared to 1/R. The interaction, then, cannot be universally 
described by an instantaneous potential and is intrinsically nonlocal in time. Thus, 
the prescription (|3|) based on Wilson loops in QCD does not yield the heavy quark 
potential in its conventional understanding. Theoretical treatment of the QQ system 
has to account for these peculiarities of QCD. Examples of the technique used here 
can be found in Ref. @ and other papers referred to there. 

In spite of all subtleties mentioned above, the heavy quark potential V defined via 
Wilson loops, Eq. (]3|), is an observable quantity (at least in the perturbative regime 
R <C 1/Aqcd); u P to an overall additive constant appearing in the renormalization of 
the straight Wilson line. In particular, it is gauge-invariant in perturbation theory. It 
is therefore tempting to use it to quantify various strong interaction effects, including 
low-energy contribution to the mass of the heavy quark. 



i b 

Figure 2: Self-energy diagram yielding the classical Coulomb shift in the mass of a 
nonrelativistic particle. 

It is well known that in classical electrodynamics the self-energy of a charged 
particle is given by the potential at the origin, ^e 2 V(0). The same holds in quantum 
electrodynamics in respect to the effect of quanta with momenta much smaller than 
the source mass (or inverse radius, for the composite particle), as is illustrated by 
the diagram Fig. 2 yielding 

Here e is charge and V em is the potential between the same-sign charges, hence it 
is positive. In QCD the perturbative diagrams for both tuq and V(q) are more 
complicated. Yet, the similar relation holds in order 0(a 2 s ) as well. This is most 
simply seen in the Coulomb gauge where all the effect at this order reduces to 
dressing the propagator of the Coulomb quanta, i.e. using the running a s (q 2 ). (An 
alternative discussion can be found in Ref. ||.) There is a general argument |l| 
showing that such a relation holds to all orders in perturbation theory: the infrared 
contribution to V(0), 




equals to minus twice the same contribution to the mass of a static color source. 

The general structure of the perturbative diagrams for heavy quark mass and 
V(0) indeed is similar. To establish the correspondence, one can cut the Q line in 
the diagram for the potential at some place and turn the new external legs around. 
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Integration over momentum q corresponds to closing the original external Q and Q 
legs, as exemplified by Fig. 3. This correspondence often works out in a nontriv- 
ial way. For example, in the potential we should discard the "reducible" diagrams 
which can be cut across only Q and Q lines, Fig. 4. a. The corresponding contribu- 
tions in SrriQ would include rainbow diagrams like in Fig. 4b. The latter, however, 
vanish in the Coulomb gauge, which is self-manifest in the coordinate representa- 
tion. The Coulomb propagator is instantaneous, S(t 1 — t 2 )/\x — y\ while the heavy 
quark propagator is retarded and includes 0{ti—t 2 ). In diagrams like Fig. 4b we 
have r 2 > T\ and therefore they vanish. The exception is one-loop diagram which 
does not vanish since in the integral over uj the large semicircle at u — > oo yields 
the finite contribution; this is equivalent to the prescription 0(0) = |. 



AAAAA 
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Figure 3: Correspondence between the diagrams for potential V(0) and for the 
heavy quark mass uiq. The dotted lines denote closing the external quark lines and 
integrating over q. 

The general argument goes as follows. Let us imagine we were able to introduce in 
some way an ensemble of gauge field configurations where the modes with momenta 
much larger than a certain scale fu, are practically absent. This field-theoretic system 
would not need regularization, and everything can be expressed in terms of the bare 
parameters, including the bare quark mass itlq\ At R — > the Wilson loop will 
approach its free value N c thus yielding V(0) = 0. This is clear on physical grounds: 
QQ form a dipole with the infinitesimal dipole moment, and its interaction with any 
soft gluon field vanishes as R goes to 0. It is important at this point that our gauge 
ensemble explicitly includes only soft modes. Otherwise, as in full QCD, the modes 
with \k\ ~ 1/R generate growing attractive potential at arbitrary small R. 
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Figure 4: Reducible QQ-on\y diagrams are not included in V(R) being iterations 
of the potential interaction. Similar diagrams in itlq vanish due to instantaneous 
nature of the Coulomb interaction. 

Next, we note that V(R) is traditionally determined up to a constant. In per- 
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turbation theory (in four dimensions) the potential is actually defined as 



U(R) = V(R) - V(oo) . (7) 

We assign U to this "standard" potential to distinguish it from V(R) which has 
a precise meaning in a finite theory. While U(R) by definition vanishes at R — > 
oo, V(oo) does not and reflects nontrivial interaction with the gluon field. In the 
momentum representation V(q) explicitly contains the term V(oo) 5 3 (q ), which is 
discarded in standard perturbative computations. Since at R— > oo the Q and Q lines 
are well separated, their interaction must vanish as 1/R, (at least, in perturbation 
theory) and the value of the Wilson loop is given by the mass renormalization of 
each static source, 

V(oo) = 2 (m^ - mg>) = 25m Q . (8) 
Therefore, we have for the "ordinary" potential 

f/(0) = / (03 V «*W = ~ 26m Q ■ ( 9 ) 

Here V reg is the regular part of V(q ) (computed in a usual way) not containing self- 
energy diagrams yielding S 3 (q). Thus, we have the stated relation between Vtr(O) 
and S m m Q . 

Relying on relation @, one can try to perturbatively define a certain running 
heavy quark mass^ which is free from the leading renormalon uncertainty ~ Aqcd 
HI. One defines 
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P Q S (H) = m p Q olc + | f V reg (q) = ml olc + I HdR V (R)[^ - fi cos ^R , 

J\q\<fi v ' JO L J 

(10) 

where V(q) is computed to a certain order in perturbation theory with no explicit 
cut-off, and the pole mass is taken to the same order in a s . The mass mo S (/i) 
is known as the "potential-subtracted" mass ||. As is clear from the preceding 
derivation, this mass would have a meaning of the rest-frame energy of the heavy 
quark, but not the mass determining the kinetic energy, which is generally different 
once a cutoff is introduced (i.e., m rather than m 2 in the notations of Refs. @, |J). 

The ansatz (pi]) is supposed to remove from the pole mass all infrared contribu- 
tions originating from the scale much below the cutoff //. As we saw, this happens 
to one and two loops. Unfortunately, this is not true starting order a 3 .. 

The problem becomes self-manifest at order where infrared-singular contribu- 

4 -*2 

tions of Fig. lb emerge in V(q) which behave like In with e being an infrared 
cutoff in the "rung" gluon momentum. On the other hand, the corresponding con- 
tributions are absent from tuq since the Coulomb exchanges are instantaneous. 

1 This idea was discussed by myself in 1996 and later, independently, advocated by M. Beneke 
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This apparent contradiction does not mean that the relation (Q) is violated. As 
explained in Ref. 0, the subtlety resides in the necessity to introduce the cut-off 
on the gluon momenta to compute Vjr. Here we consider this problem from the 
different perspective. 

Let us take a look at diagram Fig. 3a. It does contribute to V(0) computed via 
the simple-minded prescription in Eq. (|T0"D. At the same time the diagram Fig. 3c 
vanishes simply since the one-loop QQg vertex is zero for transverse gluon. The 
ansatz ( |TUD treats the diagrams in Figs. 3a and 3c differently: the integrals run 
over all momenta of the Coulomb quanta in the latter, but there is a nontrivial 
cut at momenta of order \x for potential in Fig. 3a. Our general proof actually 
states that the amplitudes of emitting the transverse gluon must vanish once the 
incoming Q and Q lines are contracted performing the integration over q. This is the 
cancellation of the diagrams in Fig. 5, which are subdiagrams to the whole potential 
V(0). Clearly, the cancellation holds only if the cuts in all diagrams are the same, 
which is not the case in the ansatz ([T0|). 

Q Q Q 
,-•-^—1 ..---5^ T- ■-■ > 
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Q Q Q 
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Figure 5: Diagrams for interaction of the transverse gluon which cancel upon 
integration over (fin the absence of hard Coulomb quanta. The cancellation, however 
does not occur for potential in Eq. (10). 

The above consideration reveals the problem of the definition of the "potential 
subtracted" mass at the technical level. In fact, there is a deeper reason behind its 
failure, which is instructive. The purpose of the additional term in Eq. (|ToD is to 
remove the infrared contributions from the pole mass. If the gluons in the diagrams 
were all soft with momenta not exceeding certain scale e, the induced potential in 
the momentum representation would vanish completely above some scale ^ > e. 
Then J\^\ <fM d 3 q V e (q) would comprise this contribution completely and correctly, 
regardless of any technical details. Why then does ansatz (|T0|) fail in practice? 

The reason lies in presence of both soft (q <^ e) and hard (q ~ /i) gluons in the 
potential, and in lack of necessary "factorization" of these scales in V(q). Let us 
single out the soft transverse gluon with momentum k ~ e. If it affected the full 
potential only at \q\ ^ e, everything would work fine. However, as exemplified by 
the H-diagrams in Fig. 1, this is not the case. Presence of hard gluons propagates 
the contribution of the soft gluon to all Fourier components of V(q). The impact 
of hard gluons on processes with soft transverse gluons is not limited to only renor- 
malizing their bare interactions, but also introduces new nontrivial 'dipole' matrix 
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elements between the heavy quark states, which, for example, do not vanish when 
the momentum k goes to zero. 

The definition of the "kinetic" running heavy quark mass based on the small 
velocity sum rules is protected against such problems by the operator product ex- 
pansion (OPE). The effect of any soft gluon on the moments of the SV structure 
functions is given by the corresponding heavy quark operators whether or not hard 
gluons are present. Likewise, the OPE bars soft physics from being transferred to 
large energy - the corresponding effects are suppressed by at least third power of 
energy || and are described by the OPE. 

The lesson one can draw from the analysis of the heavy quark potential/mass 
problem is that extreme care must be exercised in treating the infrared safe observ- 
ables for which the OPE does not apply. 

3 Heavy quark sum rules 

An important class of constraints on the hadronic parameters determining the prop- 
erties of heavy flavor hadrons follow from the heavy quark sum rules, among which 
the small velocity (SV) sum rules in the static limit rriQ — > oo play a special role. 
These sum rules are 



4 


2E l^l 2 

m 


+ 


2^ l T l/2l ' 

n 


(11) 


1 

2 


m 




2 V lr (n) l 2 

n 


(12) 


A 
2 


0Vf | T ( m )|2 
* c m| '3/2 1 
m 


+ 


Vf | r (n) l 2 

Z^ C ™I ' 1/2 1 ' 

n 


(13) 


£ = 


9 lr (m) l 2 
* tm l '3/2 1 

m 




2V f | T (n) l 2 
"1 1/2 1 ' 

n 


(14) 


3 


9V f 2 |_M|2 
^ l^i fc ml '3/2 1 
m 


+ 


2| (n)|2 
Z^ e nl r l/2l ' 

ra 


(15) 


/4 

3 


9V f 2 |_M|2 
^ fc ml '3/2 1 
m 




9 V- 2| (n)|2 
n 


(16) 


Pi 

3 


2E4l^i 2 

m 


+ 


(")|2 
Z^, t n\ ' 1/2 1 ' 

n 


(17) 


Pis 
3 


O W 3 |T-( m )|2 
^ t m 1 ' 3/2 1 

m 




oV- 3| (n)|2 
Z Z^ e nl r l/2l ' 

n 


(18) 



a sequence which, in principle, can be continued further. Here is the excitation 
energy of the k-th intermediate state ( "P-wave states" in the quark- model language), 

e k = M„ {k) - M Po , 
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while the functions t!$ and rl% describe the transition amplitudes of the ground 
state B meson to these intermediate states. We follow the notations of Ref. 0, 



(H$ /2) \A^\P q ) = -t 1/2 { Vl -v 2 ) 



2M Hr 
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(20) 



where 1/2 and 3/2 mark the quantum numbers of the light cloud in the intermediate 
states, j 77 = l/2 + and 3/2 + , respectively, and is the axial current. Furthermore, 
the slope parameter g 2 of the Isgur-Wise function is defined as 
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(21) 



Equation ( |TTD is known as the Bjorken sum rule | ID| . Super convergent sum rules 
(|12|) and ([[4]) are new [[[], [11]. Equation ([13]) was obtained by Voloshin [12]. The 
expression for kinetic expectation value /i^ is the BGSUV sum rule |L3|. The next 
one for chromomagnetic operator was derived in Ref. 



14}] , as well as Eq. 



). The 

last two sum rules are obtained along the same lines. The sum rule for the Darwin 
term p z D was first presented in Ref. 
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We have introduced the new parameter of 
the heavy quark theory E; it is the small-velocity elastic transition matrix element 
between the states with explicit spin of light degrees of freedom: for the ground-state 
vector mesons such as B* 
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(22) 

The exact magnitude of the nonperturbative hadronic parameter E is not known at 
the moment. Comparing the sum rules flT2|) , with the sum rule (p!6|) for the 
chromomagnetic expectation value //^ we expect E to be about 0.25 GeV. In the 
nonrelativistic system E is given by the product of the light mass and the orbital 
momentum, and would vanish for the ground-state mesons. 

The sum in the r.h.s. of Eq. (|H[) was previously considered in the recent paper 
||16|| . It was shown to determine one of the subleading B^D* formf actors near zero 
recoil. 

The new spin sum rules ([12]) and ([14]) are convergent and are not renormalized 
by perturbative effects. This distinguishes them from all other heavy quark sum 
rules. The unified derivation of the sum rules in the field-theoretic OPE is described 



in detail in the dedicated papers [0 , with their quantum mechanical interpretation 
elucidated. A more pedagogical derivation can be found in recent reviews ]14], |IJ. 
The new sum rules were derived in Ref. [JT^J] applying the OPE to the nonforward 
SV scattering amplitude off the heavy quark. Below we will give an alternative 
quantum-mechanical derivation. Let us, however, first discuss physics behind the 
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first sum rule (|I2p which is the relation for the total angular momentum of light 
cloud in B meson. 

At first sight sum rule fljjp which is independent of the strong dynamics looks 
surprising. In the quark models the |- and |-states are differentiated only by spin- 
orbital interaction. The latter naively can be taken arbitrarily small if the light 
quark in the meson is nonrelativistic. To resolve this apparent paradox we note 
that in the nonrelativistic case r 2 are large scaling like inverse square of the typical 
velocity of the light quark, r 2 ~ 1/Ugp. The relativistic spin-orbital effects must 
appear at the relative level ~t? S p because spin ceases to commute with momentum 
to this accuracy due to Thomas precession. The latter phenomenon lies behind the 
sum rule (1121). This connection will be elucidated below. The above relativistic 



corrections lead to the terms of order 1 in the first sum rules Eqs. (jllf), (|12|). 

To understand the connection of the sum rule (|T2|) with spin, let us consider a 
small velocity weak transition amplitude M for an elementary particle A. Let it be 
the scalar vertex, that is, mediated by the scalar current J, M—{A(v')\ J(0)\A(v)). 
We consider it in a frame moving with small velocity v. To simplify consideration, we 
can assume that the change in velocity At? is even smaller: \v'— v\ = \Av\ <C \ v\ <C 1. 
For scalar particle A the transition amplitude is described by a single formfactor, so 
that in this kinematics 



(A(v')\J(0)\A(v)} ~ const(l - a (5vv) ) . (23) 

If acceleration proceeds in the direction transverse to the velocity, this amplitude 
does not depend on the absolute velocity of the particle. 

For spin-- particle the scalar amplitude likewise is described by a single form- 
factor, however the nonrelativistic amplitudes has the new structure: 

(A(v')\J(0)\A(v)} ~ const (v?V - \ i[5vxv] -ip^aip - a (5vv ) ) (24) 

which depends on the velocity of the particle. The similar structure antisymmetric in 
v and Sv is present for any particle with nonzero angular momentum. This amplitude 
is remarkable since does not depend on the internal structure of the particle: it can 
be an elementary pointlike object, or a bound state. It depends only on its spin. 
Actual dynamics affects only the part of the amplitude symmetric in v and Av. The 
origin of the antisymmetric term is simply the transformation properties of the spin 
wavefunctions and roots to the noncommutativity of Lorentz boosts U(v ) applied 
in different directions. 

Let us illustrate this point. The state can be obtained boosting the rest- 

frame particle, U(v)\A(0)); likewise \A(v + Av)) = U(v + Av) \A(0)}. Nontrivial 
dependence of the overlap (A(v+Av)\A(v)) on the frame velocity v is nothing but 
the fact that U(v + Av) ^ U{Av) -U(v) in terms bilinear in both velocities. The 
commutator of the two infinitesimal boost operators is the rotation in the (Av, v ) 
plane; its action on a state amounts to the operator of angular momentum. Thus, 
the antisymmetric in v and At? piece of the transition amplitude for a particle 



10 



directly measures its spin.0 As is clear from the preceding discussion, the piece of 
the amplitude we are interested in comes from the phenomenon known as Thomas 
precession. 

Using this general property of the transition amplitudes, we can imagine measur- 
ing separately the heavy quark spin and the total spin of B^*> meson in the following 
gedanken experiment which would bring us close to the sum rule. We start with the 
B(*> meson at rest and consider double interaction (scattering) of the external weak 
current on the b quark, Fig. 6. In the first act the hadron is accelerated to velocity 
v, the second interaction additionally changes velocity by Av. Time t 2 —t\ between 
the two interactions is at our disposal, and we can vary it. The energy variable u 
conjugated to t 2 — ti actually measures the energy of the produced hadrons. For our 
purpose we project the final state onto B^ and, again, select only the part of the 
amplitude antisymmetric in v and Av. 




Figure 6: Double scattering process on heavy quark which measures total spin of 
B* (at large t 2 —ti) or spin of b quark (at t 2 —*ti) via Thomas precession. 

If we consider the elastic transition where the first act of interaction produces 
only moving B^*\ the amplitude is given by the total spin of flW. This chan gcs 
if we add other intermediate states. If we sum the amplitude over all b hadrons 
appearing in the intermediate state, we actually measure only the spin of the heavy 
quark. Indeed, this totally inclusive amplitude corresponds to zero time separation 
t\ = t 2 , so that the light cloud surrounding heavy quark is decoupled having no time 
to follow up the scattering. These facts are explicit in the computation of both the 
elastic contribution to the scattering amplitude, and of its OPE expansion |1|, |TTJ. 
Therefore, the contribution of the inelastic states alone directly measures spin of 
light degrees of freedom. 

Directing the interested reader to the above original papers for the formal OPE 
derivation, here we sketch the quantum mechanical way to obtain the first sum rule. 
Our starting point is the expression for the SV transition amplitudes 

(*(0)|J O (O)|SW(O)> = - Vj ^ B ^ , (25) 

and similarly for axial currents, up to spin-related factors. Here tcj = iDj is the 
momentum operator of the heavy quark. Below we also use the nonrelativistic 

2 Should we consider a vector current instead of the scalar one, there would be an additional 
term related to the Lorentz transformation of the current itself. 
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energy 7T = iD —mQ. The r.h.s. in the above relation is written in usual quantum- 
mechanical notations. In the second-quantized form 

W) W o)\#-Ho)) = _,Wg=°>lg"W°>l*-'(<>)> . (26) 

Relations for the SV transition amplitudes to the P-wave states are most easily 
obtained recalling that they are overlaps (k(v)\B(v = Q)). Then we use the general 
rule 

\H Q (v)) = \Hq(0)) -\- (v7r) \Hq(0)) + 0(v 2 ). (27) 

which nicely elucidates the meaning of the small velocity sum rules: the operator 
7r ( 7 1 ('U7r) acting on \Hq) is the generator of the boost along the direction of v . Indeed, 
to get \Hq(v)) one must find the eigenstate of the Hamiltonian with heavy quark 
moving with the momentum q = rngv. The only part which explicitly depends on 
momentum comes from the heavy quark Hamiltonian (plus higher terms in 

I/tuq). We use the relation exp (— iqx) TCq exp (iqx ) = HQ+vn+mQV 2 / 2 and drop 
the last term which is a constant; Aq obviously commutes with x. Then Eq. (p7[) 
represents the first-order perturbation theory in 67i = v7r, where — ttq plays the role 
of the unperturbed Hamiltonian Tio (further details can be found in Ref. 4 , Eq. (178) 
and Sect. VI). In the second-quantized notations the same relation takes the form 

\H Q (v)) = \H Q (0)) + J d 3 xQ^ 1 (v7r)Q(x)\H Q (0)) + 0(v 2 ) . (28) 



The first sum rule (|T2|), therefore requires computing the antisymmetric in j, k 
structure in the following sum 

^ (B*WAn) (n\n k \B*) , , 

The energies E n , Eb* are eigenvalues of the total Hamiltonian 

T~t = flight + 7~Lq (30) 



1 _o 1 



= H hght -A + - (n 2 + aB) + — ^ -(DE) + a-{Exn-nxE} +0(l/m d Q ), 

Zitlq 8m Q L 

where the I/tUq piece is the sum of the Darwin and convection current (LS) terms. 
In order to compute the sum we, following usual nonrelativistic quantum mechanics, 
use the commutation relation expressing n as the commutator of the heavy quark 
coordinate x with the Hamiltonian, but include the l/rrin terms: 

ix j = -irriQ [xj,H] + (Ej + 2i axE J . (31) 

Since each commutator with TC kills one factor E n — EB* in the denominator, the 
leading in tuq term is given by uiq ((B*\xjXk\B*) — (B*\xj\B*) (B*\xk\B*)) which 
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is symmetric over j and k. The antisymmetric part comes in the next order. To 
compute it we express, to the leading order in mq, chromoelectric field as the com- 
mutator of the heavy quark momentum and the Hamiltonian, E\ =i[7Tj, W]. Then 
we finally get 

(B-\«M(„MB-) = _ 

+ § [(B^Xk+XjTTklB*) -(^(x^-i^ixj)] (32) 
+ \ [(B*\[axn] j x k +x j [axn] k \B*)-([axn] J )(x k )-([axn} k )(x j )] . 

(The terms with separate expectation values of n and x can be discarded.) The 
basic commutator relation [xj,^] = i5j k thus yields for the antisymmetric part the 
matrix element of the heavy quark spin, | 6j k i oi . It comes from the LS term in the 
heavy quark Hamiltonian. This matrix element is given by | (e'*e k —e' k *ej), where e 
refer to the polarization vectors of B* . On the other hand, the explicit summation 
over the members of the hyperfine P-wave multiplets and their polarizations yields 
for the antisymmetric part the same spin structure (e'*e k —e' k *ej) multiplied by — t|/ 2 
and rf/2 for the §- and ^-states, respectively (see, e.g. Ref. Q, Sect. 4.1). In this way 
the general structure of the sum rule flljD is reproduced. 

A subtlety must be noted at this point which was previously discussed in the 
similar context in Ref. [17]]. In order to perform the summation over intermediate 



states, we had to include in Eq. (|32j) all possible intermediate heavy quark states \n) 
which have nonvanishing matrix elements. They can be labeled by their explicit ex- 
citation number and the total momentum (not indicated explicitly). As mentioned 
above, we consider the zero-momentum matrix elements, so that the summation 
over overall momentum is removed, cf. Eq. (p8|). This does not literally apply, how- 
ever, to the diagonal transition into the same B^ states. Although the matrix 
elements naively contain ^-function of momentum, the energy denominator in this 
case amounts to (p 2 /2i71q) 2 yielding singularity at small p. This singularity must 
be treated properly and yields a finite contribution. It must be subtracted from the 
sum, since the sum rules include only the transitions into the excited states. 

The diagonal contribution can be computed without performing actual dis- 



cretization of the problem, employing the trick used in Ref. [TjJ. Considering the 
infinitesimal momenta of B^ we can relate the matrix elements of the b quark 
momentum tt to the matrix elements of the total hadron momentum P : 

(B^(p)\n 3 \B^(0)} = ^(B^{p)\Pi\B^m ■ (33) 

In the present problem we can neglect the deviation of ratio mq/M B (*) from unity. 
We then need to evaluate the following sum over all momenta of B^: 

d 3 p (B*(0)\P j \B*W))(B*tf)\I \\B'(0)) , 



(2tt) 3 {E B *(p)-M B * 
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This, however, is nothing but the general sum rule ( |3"2"D applied to B* mesons as 
(almost) free elementary particles. (The nontrivial answer emerges due to ad hoc 
interaction assumed to render the spectrum discrete.) The antisymmetric part here 
amounts then to the expectation value of the spin of B* as a whole, and is given 
by the same polarization structure as above, but with the additional factor of 2. 
Subtracting this contribution from Eq. (^) finally yields the spin sum rule fll2[) . 

The quantum mechanical derivation of the second spin sum rule ([14]) is even 
simpler if we recall the expressions (P7|), for the moving heavy hadron state. 
On then has 

(B*(tf)|7r 3 -|B*(o))=^(B*(o)|^ 

n 

(35) 

where all quantum-mechanical states have zero momentum, and the prime indicates 
that summation does not include the diagonal transition. In writing this we have 
used the fact that the operator — 7Tq plays the role of the Hamiltonian for the states 
with zero momentum. Indeed, due to the heavy quark equation of motion tiqQ(x) = 
the relation holds for an arbitrary operator O 

(n\Qn OQ(0)\m) = (E m -E n )(n\QOQ(0)\m) . 

The transition amplitudes in the sum in Eq. (|35| ) are directly related to r's. As 
before, the contributions to the symmetric and antisymmetric in j and k structures 
there are given by ^t^, 2 +t 2 , 2 and t 2 ^ 2 — rL, respectively, which yields Eq. (|22|) with 
A and £ given by sum rules (|TB|), (|T4]). 

The spin sum rule ([12]) provides the rationale for the experimental fact that 
vector mesons B*, D* are heavier than their hyperfine pseudoscalar partners B, D. 
Indeed, if the sum rule for fi 2 G is dominated by the low-lying states then \Xq must 
be of the same sign as the constant in Eq. (0), which dictates the negative energy 
of the heavy quark spin interaction in B and positive in B* . 

The sum rules (|TTD (p~8|) obviously entail a set of exact QCD inequalities. They 
are similar to those which have been with us since the early 1980's and reflect 
the most general features of QCD (such as the vector-like nature of the quark-gluon 
interaction). The advent of the heavy quark theory paved the way to a totally new 
class of inequalities among the fundamental parameters. As with the old ones, they 
are based on the equations of motion of QCD and certain positivity properties. All 
technical details of the derivation are different, however, as well as the sphere of 
applications. The first in the series is the Bjorken inequality g 2 >\ [|H|]. We, in 



fact, have a stronger dynamical bound g 2 > | following from the sum rule (0). It 
ensures that at least some of the inelastic amplitudes must be nonzero; therefore, 
the bound state with nonzero spin of light cloud cannot be structureless (pointlike) 
regardless of bound-state dynamics. The Bjorken bound is not dynamical in this 
respect: the Isgur-Wise function — y2/ (1+vv 1 ) with Q 2 = \ is the one for the 

structureless particle. Comparison to the QCD sum rule evaluation g 2 = 0.7±0.1 
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PI suggests that the new bound can be nearly saturated. We also have the bound 
A>2£. Other bounds include p 2 > Pq and p 3 D > \p 3 LS \/2, p 3 D > —p\ s . We also 
have direct inequalities between the parameters of different dimensions, 

""-V-l' " D -8(^-l/4f ' ""-yi^TTTi)' <36) 



All these inequalities are saturated provided that only one excited state contributes. 

The old nonrelativistic quark models often do not respect the general relations 
discussed above. They typically yield small difference between the |- and |- P-wave 
states, and therefore hardly can be viewed reliable. A class of relativistic models 



was developed recently |2l| which incorporate proper Lorentz boost transformations 
and obey the exact commutation relations. Consequently, they naturally predict 
suppression of the r x /2 amplitudes compared to r 3 / 2 . Likewise, they obey the bound 
f? 2 > §. It seems appealing to reexamine the quark model predictions for processes 
with heavy flavors in the framework of the new models 

3.1 Hard QCD and normalization point dependence 



The sum rules (|lll)-(|18|) express the heavy quark parameters, including A = Mg— mj, 
p 2 and Pq, as the sum of observable quantities, products of the hadron mass differ- 
ences and the transition probabilities. The observable quantities are scale indepen- 
dent. How then, say, A = M B — m b , pi and p? G happen to be /i-dependent? 

The answer is that in quantum field theory such as QCD the sums over excited 
states are generally ultraviolet divergent due to physics at £ > Aq CD . In con- 
trast to ordinary quantum mechanics they are not saturated by a few lowest states 
with contributions fading away fast in magnitude with the excitation number. The 
contributions of hadronic states with ^> Aqcd are dual to what we calculate in 
perturbation theory using its basic objects, quarks and gluons. The latter yield the 
continuous spectrum and can be evaluated perturbatively using isolated quasifree 
heavy quarks as the initial state. The final states are heavy quarks plus a certain 
number of gluons and light quarks. It is the difference between the actual hadronic 
and quark-gluon transitions that resides at low excitation energies. 

Therefore, in order to make the sum rules meaningful, we must cut off the sums at 
some energy p which then makes the expectation values /^-dependent. The simplest 
way is merely to extend the sum only up to < p. This is the convention we 
normally use. Thus, in actuality all the sums in the relations (|TT|)-(|T8|) must include 
the condition < p, which we omitted there for the sake of simplicity, and all 
the heavy quark parameters are normalized at the scale p. The exception is the 
super convergent spin-nonsinglet sum rules (|12D and (|14D where in the perturbative 



domain p ^> Aq CD such //-dependence is power suppressed by factors a s Aq CD //x 
and can be neglected. 
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Figure 7: Perturbative diagrams determining the high-energy asymptotics of the 
heavy quark transition amplitudes and renormalization of the local operators. 



The high-energy tail of the transitions, to the order a s , is given by the quark 
diagrams in Figs. 7 with 



2E- + E 



d 3 k 

~2u 



where (u, k) is the momentum of the real gluon. The spin-singlet amplitudes are 
just a constant proportional to g s . Performing simple calculations we arrive at the 
first-order term in the evolution of, say, /i^(/i) P 



d/i 2 



4 a s 

-- + ... • 

3 7T 



(37) 



Purely perturbatively, the continuum analogs of tu^ and T3/2 are equal, and a similar 
additive renormalization of fiQ and p\ s is absent. 

The perturbatively obtained evolution equation ( |3"7D and the similar one for the 
chromomagnetic operator stating its anomalous dimension — 3a s /27r allow one to 
determine asymptotics of T\m and r 3 / 2 at e ^> A QCD , 



e m\ T 3/2 



2 I (m; 12 _^r^ 2 I ( 

fc ™! '3/2 I Z^ fc nl '1/2 
m n 



E c nl '1/2 
n 

2 u(™)|2 



ede 



8a, (e) 

9tt 
3a s (e) de 

2tt r 



E 



e t. 



(to) 1 2 



.£m<f 



ml '3/2 



£n<£ 



2 u(™)|2 
/2 I 



(38) 
• (39) 



The last bracket is |/i|.(e). Equation can be extended to higher orders in a s , 
this amounts to using the so-called dipole coupling a[ d \e) introduced in Ref. [§: 



a( d) (el 



a s (e-^ e )-3(^-g)^ + 0(a2). 



(40) 



(a s is the standard MS strong coupling). Therefore we get a number of exact 
perturbative evolution equations [|J 



dgp*) 
d/i 

dA(M) 
d/i 

d/i 2 (/i) 
d/i 



8 o^Ox) 

9 7T ; 

16 a^(/i) 



7T 



7T 



/i. 



(41) 
(42) 
(43) 
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Using Eq. ( |3Tf ) we can estimate the contribution of the high-energy states in the 
sum rules (|T2|) and ([14]): 



E 

£m<|l 


-Ml 2 
T 3/2 1 


E 


E €m 

fm<fl 


-( m )|2 
T 3/2 1 


E 

e n <n 



T- 



(n)|2 
1/2 I 

Ml 2 
1/2 I 



1 q 8 (//) AfcQf) 
4 8tt /i 2 

£ q 8 (/j) /4(/i) 

2 4tt /i 



(44) 
(45) 



they are power suppressed and presumably small in the perturbative domain. 



3.2 On the saturation of the sum rules 

The question of the saturation of the heavy quark sum rules (in particular, the 
lower ones (p~T| ) (P^6| ) ) is of primary importance for phenomenology of the heavy 
quark expansion. The sum rules at large enough normalization point fi tell us what 
the asymptotic value on the right-hand side is, and perturbation theory tells us its 
\i dependence. It is a dynamical question starting from which scale /xo this behavior 
applies. For superconvergent sum rules ( |T2"D and ( |T4D this is the question at which 
scale the sums approach the stated values with a reasonable accuracy. In order to 
sensibly apply quantitative 1/rriQ expansion, one must have mg>/i , presumably, 
mQ^>/io- While this is, probably, the case for b particles, such a hierarchy is not 
obvious a priori in charm. 

The existing numerical evaluations of A and // 2 at the scale near 1 GeV suggest 
rather large values, approximately 0.7 GeV and 0.6 GeV 2 , respectively, which impose 
rather tight constraints. These facts were often neglected under various pretexts, 
including challenging the accuracy of the QCD sum rules evaluations of yU 2 . 

Nevertheless, certain constraints following from the sum rules are tight and ro- 
bust simultaneously. Namely, the value of /z 2 ^ ~ 0.4 GeV 2 , as extracted almost di- 
rectly from B(*' and masses, has hardly been challenged. By virtue of the 
sum rules, the value of /i 2 is at least as large. Thus, regardless of the accuracy in 
evaluations of the kinetic expectation value, the question can be phrased in terms of 
the generally accepted value of (Mq. At which minimal scale /j,q the value of //^(//o) 
reaches 0.3 or 0.4 GeV 2 ? If this scale is below 1 GeV, large A and /x 2 are almost 
inevitable. If, however, GeV) is significantly below 0.4 GeV 2 , the chances for 
success in 1/tuq expansion in charm are slim. 

Similar constraints follow from the sum rule (|TJj) and, in particular, flT2"|). If these 
sum rules are saturated below 1 GeV, we expect large A and /i 2 . If the saturation 
scale is higher, perturbative treatment of the scales ~ m c only slightly exceeding 
1 GeV is not justified. 

The estimates for r 3 / 2 and T\n for the lowest P wave states group around 0.4, 

with efj 2 > e ( lj 2 ~ 400 to 500 Me V (for the review see Ref. Jy]). The old Q uark 
models yielded close values for |- and ^-states. The more recent relativistic models 
pCl predict noticeable suppression of the transitions into ^-states compared to the 
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-amplitudes. Similar absolute values were reportedly extracted for r 3 / 2 from the 



overall experimental yield of the corresponding charmed P wave states \22\. It is 



evident that such transitions amplitudes fall short of saturating the sum rules, 
n% ^ 0.2 GeV 2 , 4/2 A - °- 3 GeV > 4/2 Slight - 0-3 



(Slight denotes the sum in Eq. (|I2"D for the light cloud spin). For the spin-dependent 
sum rules, the first and the third entries, the ^-states would further decrease the 
values. In principle, the lowest states alone should not necessarily saturate the sum 
rules, even though the idea of the dominance of the lowest state contributions is very 
appealing. Let us mention that in the 't Hooft model all heavy quark sum rules are 
saturated with amazing accuracy by the first excitations p3"| |. Is such a possibility 
excluded in QCD? 

Probably, not completely. The dominance of the first excitation with e ~ 500 MeV 
(recall that one must use the asymptotic tuq — > oo values of the excitation energies 
and amplitudes) is still possible if the QCD sum rules underestimate the value of 
r 3/2-0 Experimental determinations of r's are also questionable since l/m c correc- 
tions are not accounted for there. The estimates in the 't Hooft model suggest that 
they can be very large. In the cases where they are known explicitly in QCD, the 



l/m c terms generally turn out to be very significant as well [24 



Another - and, apparently, the most natural - option is that there are new states 
with the masses around 700 MeV with similar, or even larger r^l — 0.4 to 0.5; they 
can be broad and not identified with clear-cut resonances. The |-states must be yet 
depleted up to this scale. All such states can be produced in semileptonic b decays 
and observed as populating the domain of hadronic invariant mass below or around 
3 GeV. It will be important to explore these questions in experiment. 

The branching fraction of experimentally identified narrow j '• = | P-wave states 
does not exceed a percent level. A larger fraction was reported recently for the broad 
hadronic distributions in B — > D^iriu decays [p5] , between 2% and 3.5%. Such 



an yield helps to fill the gap between the total semileptonic fraction BR S \ ~ 10.5% 
and identified exclusive channels. Since the observed D^ir mass distribution is very 
broad, these states were attributed to the j = \ resonances which are expected to 
have large decay width. If this is correct, we would rather observe significant negative 
contributions to the sum rules ([12]) and (|T6|) from the domain of e ^ 700 MeV. 
This possibility does not look natural from the theoretical viewpoint; neither it is 



supported by the predictions of the recent relativistic quark models [£0], |16| . If it is 
nevertheless realized in reality, we would have to expect quite large values of A and, 
in particular, 

We note, however, that there is an alternative interpretation of the data: the 
reported enhanced yield of broadly distributed D^n can be the manifestation of 

3 Let us note that the technology of the QCD sum rules assumes the approximate duality starting 
e = l GeV or even lower (the energies are counted from the heavy quark mass there). Therefore, 
accepting poor saturation of the exact heavy quark sum rules at this scale and simultaneously 
relying on the QCD sum rules predictions is not selfconsistent. 
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nonresonant (continuum) production; then it can belong to the |-states. If this is 
the case, one would get a consistent picture of conventional saturation of the heavy 
quark sum rules at the usual energy scale under 1 GeV, without necessity to invoke 
new paradigms like a too high onset of quark-hadron duality. 

The continuum yield is routinely assumed to be small compared to the resonant 
contribution. This, however, is more the heredity of naive quark models than the 
fact based directly on QCD. The absolute branching fraction of this yield is rather 
small, about 20% of the overall semileptonic rate, and agrees well with the general 
fact that nonresonant contributions are l/iV c -suppressed. It seems certain that such 
a possibility must be carefully explored before alarming conclusions are drawn. 

Interpreting the data mentioned above can be obscured by the potentially sig- 
nificant l/m c effects. The spin of light degrees of freedom itself is a well defined 
quantum number only as long as the limit m c — > oo is considered. While deviations 
from this academic case are possibly under control for D and D*, the situation can 
be much worse for the excited states, in particular, in the mass range we discuss. 
Here the energy of the light cloud already exceeds m c itself. The new relativistic 
quark models [ED] can gives us a sense of possible significance of such effects. 



4 BR 8 i(B) and n c 

Before concluding, I would like to mention another nonstandard possibility which 
can affect the analysis of BR s i(5) and the average number of charm quarks per B 
decay n c . As explained elsewhere, the problem of semileptonic fraction is usually 
considered in conjunction with n c - the latter constrains the total decay rate in the 
b — > ccs(d) channel and thereby allows to isolate potential uncertainties in the latter 
due to limited energy release (see, e.g., Ref. jj]], Sect. 7.1). 

The point is that there is a certain assumption in the standard analysis which 
goes beyond applying the OPE per se. Namely, the partial decay rate mediated by 
quark transition b — > ccs and computed via the OPE, is implicitly equated to the 
probability of observing hidden charm, or open charm plus anticharm states in the 
final state. Stating differently, the charmless final states and those having a single 
charm particle are presumed to be totally independent of the b — > ccs quark level 
transitions. While this is certainly the case for the single-charm channels, it is an 
additional assumption for the charmless decays; the possibility that it is violated 
was considered in Refs. [p_7| . 



Since charm quark is heavy in the scale of strong interactions, it is more than 
plausible that this assumption is accurate enough. Yet it must be realized that we 
do not derive it from the first principles of QCD. In particular, the OPE does not tell 
us this.Q The OPE predicts the total decay probability induced by the particular 
term in the weak Lagrangian, 07^(1 — 75)657^(1 — 75)0, but it cannot say if it is 

4 More precisely, I am not aware of the way to justify it from the OPE per se, without recourse 
to additional assumptions. 
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exhausted by the final states with two c quarks, or it is shared with the charmless 
final states. The OPE width must include all rescattering processes possible in the 
final state. Moreover, it is known that failure to include such contributions can even 
lead to violation of general theorem like absence of Aqcd/^q corrections to the 



decay widths ||26|| . It is worth emphasizing that we do not mean the usual effects 
from the standard Penguin operators - these properly describe the corrections to 
the total decay rates. Our point here is the effects of final states interactions which 
reshuffle the decay rate between different "apparent quark channels", but do not 
change the overall rate. 

Since charm quark is relatively heavy, the rescattering probability between cc 
and light qq (or gluonic) states must be suppressed. Yet even the relatively small 
fraction can be significant here. We can introduce the phenomenological parameter 
C to measure such a probability: 



r 6^ C c S (B -> charmless) 
C = pt5t ' ( 46 ) 



+ CCS 



where both widths are understood as generated by the the (cb)(sc) or (cb)(dc) La- 
grangians. (We assume they are normalized at the scale around m& to isolate con- 
ventional Penguin operators). By construction, C 7^ does not change the total 
nonleptonic width, but it does affect n c : 



n r — > n r 



2 C-BR OPE (B -> ccX) . (47) 



Strictly speaking, C 7^ would mandate existence of the opposite processes as well, 
say, double-charm events in the decays mediated by b — > uud transitions, beyond 
the contributions of the Penguin operators. However, the strong suppression of the 
KM disfavored channels makes such effects insignificant. 

There exists an experimental bound on such rescattering processes: BR(5 — > 
no charm) <^ 4% |28| leads to C ^ 0.2. Even being conservative, we would expect 
C smaller than that. (The estimated BR OPE (i? — ► no charm) is below 1%.) But 
the message is clear: Any value of V = BR FSI (5 — > no charm) from such final state 
rescattering would decrease n c by 2 V . Therefore, even V about 3% would have the 
dramatic impact on the allowed domain of theoretical predictions in the (BR s i,n c ) 
space. It therefore seems important to improve the direct experimental bound on 
BR(5 — * no charm) down to lto2%. 
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